Virtual annihilation contribution to orthopositronium decay rate 
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Order contribution to the orthopositronium decay 
rate due to one-photon virtual annihilation is found to be 
(5annr(2) = ( ^ ) ^ (tt^ In 0.8622(9)) Flo. 
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Positronium, the bound state of an electron and 
positron, is an excellent laboratory to test our under- 
standing of Quantum Electrodynamics of bound states. 
Although in the majority of cases the agreement be- 
tween theory and experiment is very good, the case of 
orthopositronium (o-Ps) decay into three photons is out- 
standing, since the theoretical predictions differ by about 
6 standard deviations from the most accurate experimen- 
tal result (see, however, an alternative result in [^). 
Provided that the experiment |l[ is correct, the theory 
can only be rescued if the second order correction to the 
o-Ps lifetime turns out to be ^ 250(a/7r)^rLo- 

It is however difficult to imagine how such a large 
number could appear in the perturbative calculations, 
even if the bound state is involved. This point of view 
is supported by a recent complete calculation of the 
0{a^) correction to the parapositronium (p-Ps) decay 
rate into two photons It has shown that the "natu- 
ral scale" of the gauge-invariant contributions is [several 
units] X (a/7r)^rLo- For ttiis reason, it was conjectured 
in j|] that the 0{a^) correction to the orthopositronium 
decay rate o-Ps 87 most likely is of the same order of 
magnitude. 

At first sight, the result of Ref. [Q, 

^annP^'^ = (^)' (^^ In a + 9.0074(9)) Plo , (1) 

for the gauge-invariant contribution to the o-Ps^ 87 
decay rate induced by the single-photon virtual annihi- 
lation, does not provide much support for this conjec- 



ture. In fact, the value of the non-logarithmic constant 
in Eq.(|l|) is larger by approximately one order of mag- 
nitude than the values of coefficients in gauge-invariant 
contributions to p-Ps— > 27 decay rate. 

In this note we would like to point out that the result 
for the second order correction to virtual annihilation 
contribution given in Eq.(|^) is incomplete, in that closely 
related contributions should be included as well. It turns 
out that if the missing pieces are added to Eq.(|l|), then 
the complete result for ^annP^^-' decreases and is in accord 
with the expectations advocated in Ref. 

We recall, that in bound state calculations there are 
two different types of contributions. The hard corrections 
arise as contributions of virtual photons with momenta 
k m. These contributions renormalize local operators 
in the non-relativistic Hamiltonian. For this reason they 
can be computed without any reference to the bound 
state. 

On the contrary, the soft scale contributions come from 
a typical momenta scale k ~ ma in virtual loops, and for 
this reason are sensitive to the bound state dynamics. For 
'^annP^^-' , it is casy to see that the soft scale contribution 
reads: 



4r^r(^^ = ^G(o,o), 



(2) 



where 



G(r,/) = 



\n{r)){n{r')\ 



E — En 



is the reduced Green function of the Schrodinger equation 
in the Coulomb field. 

Let us write the expansion of the Green function in a 
series over the Coulomb potential: 



G(0, 0) = Go(0, 0) + Gi(0, 0) -f G„,uiti(0, 0). 



(3) 
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The first two terms in this expansion are divergent and 
require regularization. If we use dimensional regulariza- 
tion, then the Go(0, 0) piece delivers a finite contribution, 
since it has only a power divergence. The second term 
Gi(0, 0) is logarithmically divergent. It can be easily seen 



1 



that just this term was accounted for in the calculation of 
Ref. Q , and it is precisely the term that delivers the In a 
in Eq.(|l|). However, the contributions of Go and Gmuiti 
were not calculated there. 

Both additional terms can be easily extracted from 
Ref. d. We then obtain 



A-TTfy 1 

— Go(0,0) = -a2rLo 



and 



Witir = — TT G„(o, 0) = --a^rLo, 

n=2 



(4) 



(5) 



consistent with the results of Ref.J^. 

If we now add Eqs.(|), (|) and we obtain the com- 
plete 0{a'^) correction to the o-Ps— > 87 decay rate due 
to single-photon virtual annihilation: 



<5a„nr(2) = (-) {tt^ Ina - 0.8622(9)} Flo- 



(6) 



One sees that the non-logarithmic contribution is in 
fact of order one times (a/7r)^rLOj in accord with the 
conjecture in Ref. Its value is similar in magnitude 
to the known results for other gauge- invariant 0{a^ ) cor- 
rections to the decay rate of orthopositronium [gQ. 

The only known exception from this "rule" is provided 
by the square of the 0{a) corrections to the o-Ps decay 
amplitude. This (gauge-invariant) contribution has an 
anomalously large coefficient: 28.860(2)(Q;/7r)^rLO [§)■ 
In this respect, we would like to note that there may 
be an enhancement factor due to a larger (by about a 
factor of 3) number of diagrams contributing to o-Ps de- 
cay as compared to p-Ps decay. This enhancement is 
seen already in the magnitude of the 0{a) corrections 
and translates naturally to the large value of the 0{a^) 
contribution originating from the square of the 0{a) cor- 
rections to the o-Ps decay amplitude. However, unless 
this enhancement is dramatic, it is hard to believe that 
this fact alone can explain the discrepancy between the- 
oretical and experimental results on o-Ps decay rate. 
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